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FOREWORD 

This report summarizes the results of the Phase B portion 
of Contract NAS8-27015 for the 12-month period ending 15 July 
1973. This effort was performed by the Lockheed -Huntsville 
Research & Engineering Center for NASA-Mar shall Space Flight 
Center. 

The NASA Contracting Officer's Representative (COR) for 
this study is Mr. T. C. Bannister, S&E-SSL-T. The NASA princi- 
pal investigators for this study are Mr. R.M. Poorman (M551), 
S&E-ASTN-MM; Mr. J.R. Williams (M552), S&E-PT-M; and 
Mr. E. A. Hasemeyer (M553), S&E-PT-MWM. 
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Section 1 

INTRODUCTION AND SUMMARY 

Under Contract NAS8-27015, Lockheed-Huntsville is analyzing Skylab 
Experiments M551 (Metals Melting), M552 (Exothermic Brazing), M553 
(Sphere Forming), and M566 (Ai-Cu Eutectic Growth). The primary objec- 
tive is the study of convection in the molten metals and their attendant 
solidification theory. Particular attention is given to clarifying the effects 
of reduced gravity on molten metal flow and solidification. 

Under Phase A of this contract, a ground-based study plan was prepared. 
A summary report (Ref. 1) was issued on this 10-week phase ending July 1972. 
The ground study program defined in Phase A is the subject of the present 
report. The Phase C portion of this contract will consist of a comparative 
analysis between Skylab I flight specimens and data versus those ground- 
based and KC-I35 specimens and data already processed. The Phase C seg- 
ment will continue through January 1974. The reporting of Phase B and 
Phase C results for the M566 experiments will each be delayed approximately 
three months because this experiment will be flown on Skylab IV. 

This Phase B report is concerned specifically with the Metals Melting 
Experiment, the Exothermic Brazing Experiment, and the Sphere Forming 
Experiment to be performed in the M512 facility during the Skylab I mission 
of May-June 1973. These experiments are described briefly as follows: 

M551 Metals Melting Experiment : Three sample disks, each containing 

three metal specimens of varying thicknes ses, ar e to be rotated automatically 
at a controlled speed under an electron beam gun such that an electron beam 
weld seam is produced in the metal specimen. Disk materials include 2219 
Aluminum, 321 Stainless Steel and tantalum. During the continuous weld 
portion of each disk, both full and partial penetration of the disk wili be 
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achieved by having a constant power input but a varying disk thickness. For 
each disk, the continuous weld will be followed by a dwell portion. In the dwell 
portion of the weld, the disk will remain stationary while the electron beam 
impinges on a thick segment of the disk, thus creating a large molten pool. 

The electron beam will then be shut off and the pool will be allowed to solidify. 

M552 Exothermic Brazing: A technique for joining stainless steel tubes 
will be tested in this experiment and the flow and solidification behavior of 
weightless molten braze alloys will be studied. The joining technique will 
use a solid mixture that produces heat by exothermic chemical reaction to 
braze sleeves over 0.75-inch diameter tubes, using a copper -silver -lithium 
braze alloy. A package containing four assemblies, each consisting of a tube 
with sleeve and preformed braze alloy surrounded by exothermic material, 
will be mounted in the M512 facility's vacuum chamber. The exothermic 
reactions in the four assemblies will be ignited in sequence and the whole 
package will be returned to earth for analysis. 

M553 Sphere Forming Experiment ; Twenty -eight 0.25 -inch diameter 
spherical specimens will be cast using the electron beam gun as a heat source. 
The specimens will be initially supported on two wheels by a sting. After 
melting is completed, the spheres will then be separated from their stings 
and allowed to solidify while free-floating in the vacuum chamber. Specimens 
will consist of the following materials: pure nickel, Ni-1% Ag, Ni-30% Cu, 
and Ni- 12% Sn. 

The nominal operating conditions for the electron beam gun, braze units 
and the M512 vacuum chamber are given in Appendix A. This appendix also 
contains the updated convection sensitivity analysis for all the M512 experiments 
(M55 1 , M552, M553 and M556). 

For these particular Skylab experiments, the only significant difference 
between space and earth processing will be the lack of gravity. Other environ- 
mental factors which may also differ from earth processing are the vacuum, 
radiation, electromagnetic and thermal conditions. 
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Gravity has no direct effect on grain structure or other properties of 
solidified material. These properties are determined by the crystallization 
kinetics which are controlled by short-range intermolecular forces; i.e., the 
temperature and concentration at the fluid-solid interface. Gravity has not 
been shown to have any significant direct effect on these forces, but can affect 
solidification indirectly through its direct effect on fluid motion. The three 
major indirect effects of gravity on solidification are: 

• Sedimentation 

• Bouyancy-Induced Convection, and 

• Hydrostatic Pressure. 

In addition to these three very important indirect effects, the lack of gravity 
will also afford the study of a more direct effect — the opportunity to obtain 
homogeneous nucleation. The long free-float times potentially available in 
space processing applications will allow melts to cool and nucleate without 
the deleterious effects of container walls. Wall effects are usually very 
strong in nucleation phenomena and prevent large degrees of supercooling 
from being attained in terrestrial processing. The remaining indirect 
effects are explained further below. 

Sedimentation ; This effect may be significant whenever heterogeneous 
mixtures exist in fluids, such as in monotectic, dispersed-particle, or fiber- 
reinforced composite casting. The denser of the immiscible materials will 
tend to settle unless colloidal or electrostatic attractions interfere. Also, 
in supercooled melts, segregation of freshly formed nuclei by gravity would 
affect the final grain structure. Nonmetallic inclusions, gas bubbles and 
voids, which usually exist in melts, are also distributed nonuniformly by 
g ravity . 

Convection : In terrestrial processing, gravity is the primary driving 

force for the convection of contained fluids when they are subjected to thermal 
or concentration g radients . Temperature gradients arise from external 
heating and cooling, whereas concentration g radients are usually produced 
internally (Soret effect and solute rejection at freezing interfaces which leads 
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to "constitutional super-cooling"). These two gradients can produce density 
gradients large enough to induce buoyancy-driven flow. This fluid motion 
affects the temperature and concentration profiles within the fluid. This 
subsequently alters the shape and rate of movement of the freezing interface, 
because the kinetics of freezing depend on the local temperature and concen- 
tration. The degree of mixing caused by convection may be large enough to 
change the rate of solidification from kinetic to heat transfer or diffusion 
controlled which could drastically alter the grain structure. Examples of 
this convective effect on the transitions from planar to cellular or dendritic 
growth, from columnar to equiaxed eutectic structures, and from unidirec- 
tional to colonied or banded eutectic structures are cited frequently in the 
literature. Another effect of gravity-driven convection on solidification 
processes is that the bulk fluid movement, if rapid enough, can break delicate 
dendrite arms and thereby alter final grain structure. Furthermore, inter- 
lamellar spacing in eutectic growth and dendrite arm spacing are dependent 
on cooling rate which is a strong function of convection. 

Hydrostatic Pressure : A body of fluid in a gravity field sustains a vertical 

pressure gradient as the bottom fluid must support the weight of the upper fluid. 
This pressure gradient distorts the shape of liquids on earth because the shape 
of a liquid surface is determined by the surface tension and the internal hydro- 
static pressure (and adhesion if the liquid wets a solid surface). Distorted 
drops of liquid will result in nonsymmetrical solids upon freezing. 

The following sections contain spearate discussions for each experiment 
on the effects of microgravity on molten metal flow and solidification. For 
each experiment, the following topics are discussed; the most significant 
results to date; the expected low-g variations; and specific evaluations planned 
for Phase C. In the M551 Metals Melting Experiment, the relative magnitudes 
of convection in the molten pool and their effect on grain structure and shape 
have been analyzed. In the M552 Exothermic Brazing Experiment, the effects 
of gravity on capillary flow and capillary "gap" limits have been studied. In 
the M553 Sphere Forming Experiment, the magnitude of convection has been 
determined and how the grain structure, shape and surface finish of the spherical 
specimens should be affected. 
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Visual observation of flight and ground specimens, films (M551 and M553), 
and standard post-experiment metallurgical analyses will serve as the basis 
for comparing flight test data and specimens in the Phase C study. This will 
also empirically verify the Phase B study conclusions. These observations 
will consist of identifying convective stirring from film and identifying dif- 
ferences in structure, homogeneity, and/or unidirectionality in the post- 
experiment metallurgical analyses. 
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Section 2 

M551 EXPERIMENT 

The primary objective of this experiment is to study the melting and 
solidification of metals in reduced gravity. A second objective is to evaluate 
the electron beam (EB) process as a joining and cutting technique applicable 
to assembly and repair of structures in space. The primary goal of this 
particular study is to delineate the magnitude and pattern of the molten metal 
flow for both terrestrial and space flight processing conditions. From a 
knowledge of the convective differences and the other differences between 
terrestrial and space processing, the effects on grain structure and shape of 
the weld seam and dwell are to be predicted. The "other" differences of 
space versus ground processing include differences in sedimentation and 
beading (lack of hydrostatic pressure in space) and possible differences in 
vibration, radiation and vacuum levels. These last three effects are not 
predictable to a precise degree of accuracy and are no doubt negligible in 
comparison to microgravity and therefore will not be considered important. 
Also, it was shown in Phase A studies that terrestrial versus space flight 
convective differences were negligible in the continuous weld mode (Ref. 1). 
Thus the following discussions will concentrate on the dwell mode of the 
M5 51 experiment. 

2.1 RESULTS 

A summary of the most significant results found to date are given below; 

• The following physical forces are operative during the cutting, 
joining and dwell portions of this experiment; 

Gravity/ Acceleration 

Surface Tension 

Lorentz (beam current) 

Elect restriction 

Magnetostriction 
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Electrostatic Charges 
Solidification Shrinkage 
Thermal Expansion 
V aporization 

EB Mechanical Pressure Force 

• Siginficant molten metal motion was exhibited in each of the 
specimens filmed in both the dwell and weld (disk rotating) 
modes. The films examined thus far have all been ground 
tests . 

• Both gravity and surface tension forces control molten metal 
flow during M551 ground tests. These natural convection 
driving forces arise because of severe nonis other mality during 
heating . 

• Surface tension will still provide an equal magnitude of convection 
in zero-g, although the flow pattern may be different (see Appen- 
dix B). 

• The mechanical pressure force associated with the impinging 
electron beam will not cause unstable splattering in the dwell 
mode, even in absolute zero gravity. Thus, any violent 
splattering can be attributed to surface tension or electro- 
magnetic instabilities at the weld pool surface or to degassing. 
This force is also negligible in comparison to surface tension, 
vaporization and/or hydrostatic forces during EB cutting. 

• Sedimentation will not be a factor in this experiment, because 
vigorous stirring will occur both in terrestrial and Skylab 
operations . 


The preceding conclusions were drawn from examination of color motion 
picture film taken during ground tests and from detailed examination of the 
descriptive differential equations which govern EB melting and solidification. 
The equations were examined using dimensional analysis (Refs. 2 and 5 ) and 
the nominal operating and boundary conditions for this experiment and the 
M512 chamber (Appendix A). Detailed derivations and results are given in 
Appendix B . 
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It should also be mentioned that surface tension, gravity and the electro- 
magnetic forces are the most significant of all the preceding physical forces 
for convective motion. On Earth, gravity and surface tension art' equal and 
dominant; whereas, on Skylab, surface tension will be dominant, followed by 
electromagnetic, then by gravity. Furthermore, there are two primary 
effects of surface tension on melt behavior; 

• The mean effect of uniform surface tension results in normal 
stresses at free surfaces which determine the shape of fluids. 

• Nonuniform surface tension leads to tangential stresses at free 
surfaces which can induce fluid motion. 

2.2 EXPECTED LOW-G VARIATIONS 

Based on the preceding results and conclusions, the following predictions 
are made regarding the variations expected in the disk welds processed in the 
microgravity conditions aboard Skylab; 

• Slightly different grain structure and/or dendrite arm spacing 
may result due to different melt flow patterns (see Appendix B). 

• Lack of hydrostatic pressure and dominance of surface tension 
will allow beading and dwell pool shape to be different (more 
spherical) and may cause deeper EB penetration. 

The latter will occur because surface tension acts to pull the melt toward 
cooler surfaces which is away from the electron beam in this case. Assuming 
that the net gravitational field is always aligned with the beam in Skylab 
processing, the situations illustrated in Fig. 2-1 will result regarding bulk 
melt behavior of the dwell pool. As seen in Fig. 2-la, surface tension domi- 
nance will pull the melt away from the beam, resulting in less resistance and 
deeper beam penetration than in terrestrial processing when healing from 
above. When heating from the side, however, terrestrial processing allows 
the melt to run out of the pool which will probably lead to even less resistance 
and more penetration than in space processing. Furthermore, the cut or 
seam may be asymmetric because of the uneven distribution of the melt in 
the terrestrial heating -from-the- side configuration. 
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2.3 PHASE C EVALUATIONS 

Specific evaluations to be performed for this experiment in the Phase C 
portion of the study will include: 

• The flight film will be examined and compared with ground test 
films to investigate differences in magnitude (velocity) of con- 
vection and flow patterns. 

• The microstructure and grain size (to be obtained from other 
investigators) will be examined and compared with ground test 
results to observe if different melt shapes and convection 
patterns had any significant effects. 

In addition, visual observations of flight and ground samples will be made to 
determine if beading and penetration changed significantly in space processing. 

2.4 CONCLUSIONS AND SUMMARY 

The primary conclusion of the preceding study is that significant stirring 
of the melt will occur both in terrestrial and space processing for the M551 
experiment. The magnitude of this convective fluid flow will be the same in 
both environments. The flow patterns, however, will be different because 
surface tension alone will be the convective driving force for Skylab proces- 
sing whereas gravity and surface tension forces will be equally coupled driving 
forces on Earth. As a result of surface tension dominating pool shape, 
quicker, cleaner, narrower and more symmetrical cuts and penetrations 
should be attainable in EB melting in space. 

It is recommended that the M551 experiment be repeated in the space 
environment essentially as is except that high-speed photography should be 
added. This will allow important hydrodynamic phenomena in cutting, bead- 
ing, and the dwell to be studied. The current 24 ft/sec filming rate of the 
flight test will miss much of the hydrodynamic behavior of the melt. An 
accurate description of the hydrodynamics involved may settle debates which 
currently exist over how the beam penetrates and what causes rippling. 

Further recommendations include the following additions to M551 experi- 
ment objectives: 
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Study surface tension dominated fluid flow in the absence of 
gravity. In terre strial studies, gravity forces usually pre- 
dominate and cloud all other effects. 

Study the fluid dynamics involved in EB melting. Understand- 
ing of the stirring and heat flow associated with convective 
motions should be a valuable contribution to EB welding re- 
gardless of whether terrestrial or space applications are 
being considered. 


11 
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Section 3 

M552 EXPERIMENT 

The primary objectives of this experiment are to study and evaluate the 
capillary flow of molten braze material, to evaluate a tube joining technique 
for the assembly and repair of hardware in space, and to demonstrate the 
feasibility of exothermic reaction in space. The primary goal of this par- 
ticular study is to predict the effects of gravity on capillary flow and capillary 
"gap" limits. Accounting for these differences in capillary flow and other 
space-ver sus- Earth processing differences (e.g., sedimentation, lack of 
hydrostatic pressure, vacuum, etc.), changes in braze solidification are to 
be estimated. 


3.1 RESULTS 

A summary of the most significant results found to date are: 

• Spreading time will be up to 50% shorter in gravity fields antici- 
pated for Sky lab than in ground tests. 

• Turbulent or oscillatory laminar flow will occur for some gap 
widths in M552 Skylab processing, whereas only steady laminar 
flow would occur on Earth. 

• Marangoni bubble migration will probably be important for M552 
Skylab conditions, whereas buoyancy will control bubble location 
under terrestrial conditions. The Marangoni effect is a surface 
tension force which can drive fluid motion. It can arise from 
temperature or concentration gradients. In the present case, 
only temperature gradients are considered. 

• The extent of capillary rise or spreading will be increased in 
zero gravity. This will be an especially significant result in 
cases where the capillary driving pressure has the hydrostatic 
pressure of the liquid column opposing it during terrestrial 
operations . 

• The density ratio of silver to copper is p^ ,/p =1.2, which 

means that silver may have a tendency to s^ttle^lo the bottom 
of the joint on Earth, but would remain mixed on Skylab. 

12 
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The detailed calculations and data upon which the preceding results were 
derived are given in Appendix C. 

It is very important that the spreading time and fluid flow conditions 
be known in brazing. At first glance the flow analysis suggests that the 
time required to fill the joint is insignificant compared with the overall time 
of the brazing operation. This is not the case, however, when one considers 
rapidly occurring interactions such as alloying between the base and filler 
metal, differential rates of vaporization between the braze alloy constituents, 
and electrocapillarity (layers of ions alter the surface free energy of the 
filler). 

These and other detrimental interactions, along with the importance of 
a knowledge of fluid flow conditions, are explored in detail by Milner (Ref. 4). 
Many of the important interactions with capillary flow in brazing operations 
are due to the presence of the fluxing agent. Fluxes remove the initial oxide 
film but may not entirely prevent further reoxidations. The use of a flux 
may also modify wetting and filling of a joint by: 

• Electrochemical reactions between the base and filler metal 
with the flux acting as an electrolyte. 

• Electrocapillarity in which layers of ions change the surface 
free energy. 

• Substitutional chemical reactions in which an element of the 
flux is deposited on the surface of the base metal. 

Alloying between the braze and base metal can radically affect joint filling by 
changing the surface tension, viscosity and melting point of the filler metal. 
Surface roughness can have an important effect as it may reduce a contact 
angle of up to 80 degrees down to zero. Surface roughness is also related 
by some to the well-known dynamic capillary effect known as contact angle 
hysteresis (Refs. 5 and 6). The quicker the molten braze spreads, the less 
time there will be for the preceding interactions to significantly alter 
spreading. 
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3.2 EXPECTED LOW-G VARIATIONS 


Using the preceding results and conclusions as a foundation, the follow- 
ing P rec *i c tions can be made regarding the variations expected in the brazed 
joints processed in the microgravity environment aboard Skylab; 

• The extent, rate and uniformity of spreading should be increased 
in low gravity which should lead to better joints. If turbulence 
occurs in wider gaps, this improvement may not be attained. 

• The lack of hydrostatic pressure should allow a more symmetrical 
interface at the free surface of the filler and permit capillary 
spreading in wider gaps . 

• The location of bubbles or voids should be different between 
Earth and Skylab processing and will depend on container 
orientation to gravity direction and the magnitude of effective 
gravity during the brazing operation. 


The first item above is based on the premise that faster, laminar flow 
rates will allow less time for deleterious interactions to occur in the spreading 
molten braze. 


The last two variations are directly dependent on the exact magnitude 
and direction of the gravity vector which will be operating during Skylab 
processing of M552. The magnitude of gravity aboard Skylab during the 
M552 experiment will be approximately 10”^ g^, (gg = 9.8 m/sec), but the 
gravity direction is not yet known. Thus, precise predictions on bubble or 
void location cannot be given at the present time. 

3.3 PHASE C EVALUATIONS 

Specific evaluations to be performed for this experiment in the Phase C 
portion of the study will include: 

• The extent and uniformity of spreading will be examined using 
data supplied by other investigators. 

• The location of voids will be evaluated taking into account 
solidification shrinkage, buoyancy, Marangoni effect and the 
presence of the tube slot. 

14 
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• The microstructure data obtained by others will be analyzed, 
especially near the free surface where surface tension con- 
vection might occur even in low gravity. 

In addition, visual observations of sectioned flight and ground samples may 
provide direct information of the extent of spreading and voidage changes 
which occurred in space processing. 

3.4 CONCLUSIONS AND SUMMARY 

The primary conclusion to be drawn from the preceding discussions is 
that better quality joints may result due to the increased extent, rate and 
uniformity of spreading attainable in the microgravity environment aboard 
Skylab. If turbulence occurs in the wide gap specimens, this improvement 
may not be seen for these specimens. Also, based on the aforementioned 
results and conclusions, a new experiment involving the brazing of extra- 
wide gaps and several new configurations of variable size gaps (e.g., 
eccentric annuli at acute angles of departure) is suggested. 
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Section 4 

M553 EXPERIMENT 

The basic objective of this experiment is to determine the effects of 
reduced gravity on fundamental solidification phenomena. In particular, 
the following two effects should be prominent: 

• The density difference between components 
(e.g., nickel and tin) will not lead to sedimen- 
tation as might occur in terrestrial opera- 
tion. 

• The specimens may be cast from the melt 
without using containers or molds (homogeneous 
nucleation). 

The primary goal of this particular study is to delineate the magnitude 
and pattern of the molten metal flow for both terrestrial and space flight pro- 
cessing conditions. From a knowledge of the convective flow differences and 
the other differences between terrestrial and space processing, solidification 
theory will be utilized to predict how the grain structure, shape, internal 
voidage and surface finish of the spherical specimens should be affected. 

The "other" differences of space versus ground processing include changes 
in sedimentation and sphericity (lack of hydrostatic pressure and dominance 
of capillary forces in space) and possible differences in vibration, radia- 
tion, electromagnetic, and vacuum levels. The last four effects are not 
predictable to a precise degree of accuracy and are no doubt negligible in 
comparison to microgravity effects and will therefore not be considered im- 
portant. 
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4.1 RESULTS 

A summary of the most significant results found to date are given below. 


9 The following physical forces are operative on the 
M553 specimens during EB melting and solidification: 

• Gravity/Acceleration 

• Surface Tension 

• Lorentz (beam current) 

• Electrostriction 

• Magnetostriction 

• Electrostatic Charges 

• Solidification Shrinkage 

• Thermal Expansion 

• Vaporization 

• EB Mechanical Pressure Force 


• Vigorous convective stirring occurs in the molten 
metal for each of the M553 sample materials. This 
fluid motion was exhibited on high speed motion pic- 
ture film taken for both ground tests and low gravity 
KC-135 aircraft tests. 

9 Both gravity and surface tension forces control molten 
metal flow during M553 ground tests. These natural 
convection driving forces arise because of severe non- 
isothermality during heating. 

9 Surface tension driven convection will occur in Skylab 
electron beam melting. In the near absence of gravity, 
the surface tension forces will provide an equivalent 
amount of convection aboard Skylab, but the flow pattern 
may be different from those of operations on Earth. 

• Velocities of 20 cm/ sec magnitude will be attained 0.1 
sec after melting begins in the M5 53 nickel specimen 
(both for Skylab and ground tests). Fluid flow will 
decay 60 seconds after melting begins. This means 
that some degree of flow will exist upon freezing even 
for the free-floating M553 spheres. 


17 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306697 


• The mechanical pressure force associated with the 
impinging electron beam is not the cause of violent 
droplet breakup exhibited in some of the earlier KC- 
135 aircraft tests. Thus, any unstable splattering 
of the melt during heatup must be attributed to sur- 
face tension or electromagnetic instabilities at the 
melt surface or to degassing. 

• Sedimentation will not be a factor in this experiment, 
because vigorous stirring will occur both in terres- 
trial and Skylab operations. 


The preceding conclusions were drawn from examination of color motion 
picture film taken during ground and low-gravity aircraft flight tests and from 
detailed examination of the descriptive differential equations which govern EB 
melting and solidification. The equations were examined using dimensional 
analysis (Refs. 2 and 3) and the nominal operating and boundary conditions for 
this experiment and the M512 chamber (Appendix A). Detailed derivations and 
results are given in Appendix D. 


It should also be mentioned that surface tension, gravity and the electro- 
magnetic forces are the most significant of all the preceding physical forces 
for convection. On Earth, gravity and surface tension are equal and dominant; 
whereas, on Skylab, surface tension will be dominant, followed by electromag- 
netic, then by gravity. Furthermore, there are two primary effects of surface 
tension on melt behavior: \ 

• The mean effect of uniform surface tension results in 
normal stresses at free surfaces which determines 
the shape of fluids. 

• Nonuniform surface tension leads to tangential stresses 
at free surfaces which can induce fluid motion. 

4.2 EXPECTED LOW-G VARIATIONS 

Based on the previous results and conclusions, the following predictions 
can be made regarding the variations expected in the specimens processed in 
the microgravity environment of Skylab: 
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• The primary differences expected in Skylab samples 
will be due to containerless casting; i:e., there will 
be no wall effects upon nucleation and solidification. 

Thus dendritic and grain structure should be finer 
because of the large degree of supercooling which 
can be obtained in homogeneous nucleation. Contain- 
erless, free suspension will also provide more per- 
fectly spherical castings. 

41 Solidification shrinkage may induce imperfections in 
the sphericity of the castings. Furthermore, shrink- 
age may cause a number of voids or one large shrink- 
age cavity in the interior of each casting. The latter 
may occur symmetrically and form the basis for "hollow 
ball bearings." 

• Finally, the residual fluid motion existing at the instant 
of freezing may introduce some inhomogeneity in the 
microstructure. 


4.3 PHASE C EVALUATIONS 

Specific evaluations to be performed for this experiment in the Phase C 
portion of the study are given below. 


• The flight film will be examined and compared with 
KC-135 aircraft and ground test films to investigate 
differences in magnitude (velocity) of convection and 
flow patterns. 

• The microstructure and grain size (to be obtained 
from other investigators) will be examined for evi- 
dences of convective effects. 


In addition, visual observations of flight and ground specimens will be performed 
to determine if sphericity and surface finish improved markedly in space processing. 


4.4 CONCLUSIONS AND SUMMARY 


The primary conclusion of the preceding discussions is that significant flow 
and stirring in the molten spheres (during and after EB melting) will occur in both 
terrestrial and space processing. The magnitude of this convective fluid motion 
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will be the same in both environments. The flow patterns, however, will be 
different because surface tension alone will be the convective driving force 
for Skylab processing whereas gravity and surface tension forces will be 
equally coupled driving forces on Earth. As a result of surface tension dom- 
inating specimen shape, more perfectly spherical and smoother spheroids 
should be obtainable. This experiment may also engender the production of 
thick-walled hollow spheres which would have ready application as hollow 
ball bearings. Furthermore, the unique microstructures attainable through 
containerless nucleation and solidification may yield castings with unusual 
structural (strength) or electromagnetic properties. 

It should be noted that a repeat of this experiment in the microgravity 
environment of space as is, except for longer, higher speed camera coverage, 
would yield much more hydrodynamic information. This information is im- 
portant because hydrodynamics controls the release of the specimens and 
affects the solidification and mic rostructural properties. The release mech- 
anism is of utmost importance because it will control the duration of free 
float time in the absence of acoustic or electromagnetic levitation. It is 
recommended that a new experiment be performed with no active mechanical 
forces (springs) for release, i.e., allow only the vaporization and electro- 
magnetic forces to initiate release. The slight forces could be counteracted 
by the small residual gravity field to yield very slow release acceleration 
and extremely long free-float times. 

Further recommendations include the following additions to M553 exper- 
iment objectives. 

• Study surface tension driven convection (due to non- 
uniform temperature) in free floating, molten spheroids. 

Similar types of fluid flow probably exist in stars such 
as the Sun, and this information should be of great 
interest to astronomers. 

• Study the hydrodynamics and kinematics of deployment 
in zero-g. This step is critical for obtaining long 
periods of free float which will be essential to future 
space processing applications. 
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• Study the fluid dynamics involved in electron beam 
melting. Understanding of the stirring and heat 
flow associated with convective motions should be 
a valuable contribution to EB melting and solidifi- 
cation regardless of whether terrestrial or space 
applications are being considered. 
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Appendix A 

In this appendix a compilation is presented of pertinent environmental 
conditions for the M512 experiments to be performed aboard the first Skylab 
mission. Also included is a summary of convection sensitivity calculations 
for the M 51 2 experiments. 

The M512 facility will contain several experiments (M551, M552, M553, 
M566) for the investigation of melting and solidification in a space environ- 
ment. Convection in the molten fluid may be a controlling parameter in these 
processes. In order to evaluate convection and other phenomena, the pertinent 
environmental conditions of each experiment must be known. This appendix 
contains a compilation of the factors that are important for assessing convection. 

The environmental conditions expected for each experiment are shown in 
Table A-l. Data enclosed within parentheses indicate that they are estimates 
and need further study. Also some categories have no data which indicates 
that these factors could not be estimated at the present time. 

Dimensionless ratios which give an indication of convection sensitivity 
for each M512 experiment are shown in Table A-2. These numbers were 
evaluated by estimating pertinent temperature gradients from documented 
ground-base studies, Brown Engineering thermal analyzer studies and Lock- 
heed thermal studies and utilizing the data in Ref. A-l, The table contains 
Rayleigh numbers (Ra), Marangoni numbers (Ma), thermal diffusivities (a), 
Prandtl numbers (Pr), and Stefan numbers (St) which have been obtained for 
the molten metals in the M551, M552, M353 and M566 experiments. Critical 
Rayleigh (Ra c ) and Marangoni (Ma c ) numbers for these experiments were 
estimated from the literature (Refs. A-2, A-3, and A-4). When the Rayleigh 
number exceeds its critical value, gravity-induced convection will occur and 
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Table A-l 

ENVIRONMENTAL CONDITIONS FOR M512 


Condition 


Date/Time of Experiment 

2 

Gravity Level (cm/sec ) 

2 

Gravity Gradient (l/sec ) 

Initial Chamber Pres, (atm) 

Initial Atmos. Composition 
in Chain be r (wt %) 

Initial Chamber Temp. (°C) 

Electron Beam 

• Wattage 

• Amperage 

• F ocus 

• Duration 

• X-ray level 

Magnetic Field Level 
Vibration Levels 

• Amplitude (cm) 

• Frequency (Hz) 

Electrical Fluctuations 
Wheel Rotation 

• Rate (rpm) 

• Radius (cm) 

• Weld Radius (cm) 

Cabin Conditions 

• Temperature (°C) 

• Pressure (atm) 

• Composition (wt %) 

Space Atmosphere 

• Temperature (°C) 

• Pressure (torr) 

• Composition 

Chamber Size 

• Volume: 

• Diameter: 

Size of Valve to Space: 


M551 


M552 


M553 


M554 


1.6 kW 
80 mA 
.03 - .09 in. 
15-30 sec 
(> 20 jLirad/ sec) 


2.5 + 0.1 
8.25 
6.04 


r 5/8/73 
(.001 - . 1 ) 

I 

(0.34) 

(O 2 (70), N 2 (30) 

L (10 - 32) ^ 


(1.2 kW) 
50-60 mA 
.25 in. 

3-5 sec/sphere 
(>20 jurad/ sec) 


20 - 25 
0.34 

O 2 (70), N 2 (30) 


I 950 - 1160 

L- o l 10“ 6 - 10" 8 J 

10 e/cc, 10° neut. particles/cc, MW = 20-26 


1.5 cu ft 

16.25 in. 

4 in. vent valve (40 in. long tube to space) 
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significantly affect heat transfer. Similarly, if the critical Marangoni number 
is exceeded, surface tension-driven convection will occur. The Stefan number 
is useful in analyzing radiative cooling (M551 and M553) and the thermal 
diffusivity and Prandtl number are useful in studying thermal profiles and fluid 
motion. Table A-3 represents a quick comparison of the modes of convection 
expected in M512 experiments on Earth and on Skylab. 

/ 

It must be emphasized that the critical conditions only rigorously apply 
when gravity is aligned with and opposed to the temperature gradient. It is 
also assumed that uniform heating exists (linear temperature gradient) and 
that temperature gradients only occur in one direction. With these restrictions 
in mind, the conclusions in Table A-3 give an accurate indication of the con- 
vection modes which should be expected for these experiments. 

No values for Ra and Ma are shown for tantalum because a reliable 
estimate of tantalum's thermal conductivity in the molten state has not yet 
been found. Similarly the critical Marangoni numbers for M566 and M552 
are not given because Ma c has not been estimated in the literature for such 
large aspect ratio configurations. Values of the Stefan number for these two 
experiments are not given because melt radiation will not be a factor. 

The Rayleigh number, Ra, is the weighted ratio of the buoyancy driving 
force to the viscous forces which tend to resist fluid flow and is defined by, 

Ra = 

VCL 

where g is the gravity acceleration, d the depth of the fluid layer, J3 the 
thermal volumetric expansion, 4T a temperature difference in the fluid, 

V the kinematic viscosity, and a the thermal diffusivity. A corresponding 
number of surface tension force to viscous force is given by the Marangoni 
number: 
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Table A- 3 

GRAVITY AND SURFACE TENSION CONVECTION IN M512 EXPERIMENTS 




Gravity - 

Driven 

Surface Tension- 
Driven 

Experiment 

Material 

lg 

o 

i 

OP 

M55 1 

2219 A l 
(Rotating Disk) 

No 

No 

No 


321 SS 
(Dwell) 

Yes 

No 

Yes 

M552 

Ag-Cu Braze Alloy 

No 

No 

No 

M553 

Nickel 

Yes 

No 

Yes 

M554 

Af-Cu 

Yes 

No 

No 
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Appendix B 

This appendix contains detailed results of calculations for the M553 ex- 
periment which were performed during the Phase B portion of this study. The 
first section contains an analysis of the driving fore es, magnitude and pattern 
of natural convective fluid flows which occur in the Sphere Forming Experi- 
ment. This section is followed by a discourse on the possibility of molten 
pool splattering due to the momentum of the impinging electron beam. 

B. 1 CONVECTION ANALYSIS 
B. 1 . 1 Summary 

Results to date of the convection analysis indicate that: (1) surface 

tension is the main driving force for fluid motion in microgravity electron 
beam welding; (2) gravity forces also predominate in ground-based welding; 

(3) vapor pressure is the primary cause of surface deflection (cutting action); 

(4) high-speed photography (2000 frames per second) of KC-135 or ground- 
based tests are required; and (5) fluid dynamics undoubtedly control the 
mixing, microstructure, and shape of the weld pool in the electron beam 
weldings. 

B.1.2 Dimensional Analysis 

As in the other M512 experiments, the fluid dynamics of the M551 ex- 
periment is the most important factor in determining quality of the final pro- 
duct (in this case, the solidified pool, seam and cut produced on various areas 
of thin metal disks). Flow patterns in the molten material are important in 
these experiments, because all of the materials have low entropies of fusion 
(Ref. B-l). Thus, their solidification (microstructure) is controlled by the 
rate of heat transfer removal (Ref. B-2), which changes with the fluid flow 
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(Refs. B-3 and B-4). The flow will be especially important in the dwell mode, 
since a relatively large pool of melt will be created. The degree of flow will 
also determine the amount of mixing attained. If no or little flow were present, 
all heavier components would segregate to the bottom of weld zones on earth, 
but to a much lesser extent in microgravity environments. Fluid flow can also 
affect the shape of the weld pool (Ref. B-7). 


Application of dimensional analysis (Refs. B-5 and B-6) to the governing 
equations for EB welding, coupled with ground-based and KC-135 experiments, 
should enable prediction of the entent of reduction or increase of motion in the 
weld pool and/or the change in flow pattern in electron beam welding in space. 
Possible physical forces which could induce fluid flow in the M551 experiment, 
and their causes, include: 


• Effective Gravity Force: Resultant force on weld 

specimen due to earth* s gravity and centrifugal and 
coriolis forces of orbiting spacecraft. 

• Lorentz Force: Electromagnetic forces induced by 

passage of the electron beam current through the 
specimen. 

• Elect restriction: Stresses induced when electrical 

permittivity changes with density. 

• Magnetostriction: Stresses induced when permeability 

changes with density. 

• Electrostatic Force: Caused by presence of excess 
electrical charge (due to beam current and/or thermionic 
emission). 

• Surface Tension : Tangential stresses at vapor-liquid 

or liquid-liquid interfaces can be induced if surface 
tension depends on temperature and/or concentration. 
Surface tension will also cause pressure gradients 
across curved interfaces. 

• Density Differences Accompanying Phase Changes 

• Beam Force : Impinging electrons give up their 

momentum. 
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• Thermal Expansion: Dilation and compression of 
fluids whose density changes appreciably with temper- 
ature can induce fluid flow. 

• Vibration: Uncontrolled movement due to engine opera- 

tion, astronaut motion, particle impacts, etc. 

• Centrifugal and Coriolis: Generated by disk rotation. 

• Vapor Pressure: Evaporating molecules impact 

momentum which leads to normal stresses at vapor- 
liquid interface. 

• Inertia Forces: Tend to sustain induced motions. 

• Viscous Forces: Tend to resist driving forces. 


The preceding forces, which could influence fulid flow and solidification, 
appear explicitly in the conservation equations which apply to formation of a 
molten pool by electron beam heating. These equations are given in detail in 
Appendix D. 

A formal method of determining the controlling physical forces affecting 
fluid flow and solidification in electron beam welding was introduced in Appendix 
D, Section D. 1„2„ The controlling physical forces are determined by non- 
dimensioning the governing differential equations and performing an order- of- 
magnitude comparison on the various dimensionless groups which result. The 
key to successful analysis is in choosing the proper reference values; i.e., 
since no freestream velocity exists, which forces do we equate to estimate a 
"characteristics" or "typical" velocity. Choosing the proper characteristic 
velocity is very important, since the reference time, temperature, etc., usually 
depend on this velocity. The results arrived at in Appendix D, Section D.1.2 
are also valid for the M551 experiment. 

The equation, in dimensionless form, governing electron beam melting 
is Eq. D-l shown in Appendix D. Values of the pertinent dimensionless groups 
for each of the M551 materials are given in Table B-l. Examining Table B-l 
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Table B-l 

M551 DIMENSIONAL ANALYSIS 





1 

❖ 

N Bo 

N Bo + 

Material 

N Oh 

N q 

Bo 

N Oh 

N Oh 

N Oh 

2219 Aluminum 

5.8 x 10" 4 

3.09 

1713 

5287 

0.53 

321 Stainless Steel 

6.3 x 10' 4 

1.56 

1581 

2461 

0.25 

Tantalum 

3.4 x 10' 4 

7.19 

2942 

21150 

2.11 


With Earth gravity (9-8 m/s ec ) 
^With 10"^ Earth gravity 
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in conjunction with Eq. (D.l), an order- of-magnitude analysis indicates that 
surface tension driven convection will occur both in ground tests and for Skylab 
conditions because Nq^CCI. Furthermore, gravity driven convection will exist 
on ground tests, but will be negligible in the reduced gravity of Skylab. Thus 
different forces will control convection on earth versus Skylab. The preceding 
analysis also indicates that electromagnetic or Lorentz forces will be negligible 
with regard to causing fluid motion. See Appendix D, Section D. 1.2 for further 
details. 

A review of the literature on electric arc welding, which is somewhat 
similar to electron beam welding, has generated the following facts which can 
be compared to the preceding conclusions. In a study of motion in weld pools 
in arc welding, Woods and Milner (Ref. B-7) conclude that Lorentz forces are 
the primary cause of motion. An additional secondary cause is the momentum 
imparted by the impinging arc. Surface tension forces were not considered. 
Kotecki et al. (Ref. B-8) showed that surface tension and momentum forces 
of the impinging arc controlled ripple formation in gas tungsten arc welds. 
Surface tension controls when the arc is shut off. Brimacombe and Weinberg 
(Ref. B-9) also conclude that surface tension is the driving force for fluid motion 
once the impinging jet is removed. 

One direct result of formulating the governing equations is that surface 
deflection (cutting action) is primarily caused by vapor pressure, but that the 
beam force is also appreciable. This was determined from the vapor-liquid 
force equation of Section D.1.2, which indicates that, 


Vapor Pressure 
Beam Pressure 


This agrees with earlier studies (Ref. B-10). 
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B.1.3 Flow Patterns 

The aforementioned difference in controlling forces between terrestrial 
and space processing for the M551 experiment should result in altered flow 
patterns between these two conditions. From a knowledge of flow patterns 
generated by surface tension versus buoyancy controlled convection on Earth 
(Ref. B-ll), the situation illustrated in Fig. B-l would no doubt exist during 
melting on Earth, while that shown in Fig. B-2 would prevail in space. These 
two figures represent cross-sectional views through the molten dwell pool. 

The pattern shown in Fig. B-2 is abased on the fact that surface tension is 
stronger in cooler regions; therefore the coolest portion of the specimen will 
exert the strongest "pull" on the surface. In both cases, the motion is termed 
cellular (Ref. B-ll). 

These predicted flow patterns have been verified by M551 ground tests. 
Movies taken during ground-based tests of the M551 experiment for all three 
materials (2219 aluminum, stainless steel and tantalum) were studied during 
the week of 16 March 1973. The quality of the film and the resolution of the 
fluid flow in the weld pool were excellent in several of the tests. Four distinct 
Benard convection cells were observed in the weld pool of a stainless steel 
dwell. Significant molten metal motion was exhibited in each of the specimens 
filmed in both the dwell and weld (disk rotating) modes. 

3.2 WELD POOL SPLATTERING 

In both the M551 and M553 experiments, the heat released to the metal 
specimens by the impinging electron beam enables a molten pool to form. A 
question arises as to the stability of these molten liquids in low gravity. At 
least one of the KC-135 M553 specimens was seen to break up violently into 
many smaller liquid spheroids upon complete melting. The electron beam 
was still hitting the specimen during break up. A similar instability might 
develop in the molten puddle formed during the dwell mode of the M551 exper- 
iment, wherein the liquid might not adhere and separate from the solid disk 
(violent splattering). 
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Fig. B-l - Convective Pattern in M551 Dwell Pool 
Melting on Earth 



Fig. B-2 - Convective Pattern in M551 Dwell Pool 
During Melting in Space 
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There are at least three different mechanisms which might explain the 
instability observed in the M553 KC-135 specimens. Upon melting, the spec- 
imen may have experienced violent degassing which could have led to droplet 
breakup. Another mechanism might be electrohydrodynamic instability caused 
by interactions between the electromagnetic forces of the electron beam and 
the fluid flow field set up in the molten metal by both thermal gradients, Lorentz 
and other forces {Ref. B-12). Lastly, the momentum force associated with the 
impinging electron beam might have set up unstable surface oscillations on the 
molten metal. 

The latter instability mechanism has been treated recently by Berghmans 
(Ref. B-13). He performed a theoretical study of fluid interface stability with 
special attention being given to the role of surface tension. It was motivated 
by its possible application to the splattering of molten metal as observed during 
electric arc welding. 

Berghmans' study concluded that the weld pool interface would be stable 
if the following condition was met, 

We 2 < 1.04 + 3.3 Bo 2 (B.l) 

where 

We = Weber number 

Bo = Bond number 

The analysis is only rigorous if inertia and viscous effects are small compared 
to surface tension effects. Inertia effects are not negligible in the M551 and 
M553 experiments, but the results of Berghmans' study should give a reasonable 
approximation. 

The criterion expressed by Eq. (B.l) were applied to each of the materials 
in both the M551 (2219 aluminum, stainless steel, and tantalum) and the M553 
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(pure nickel, Ni-Cu, Ni-Ag, and Ni-Sn) and for gravity levels between ground 
tests (g-g = 980 cm/ sec ) and those of Skylab (g = 10 g^,). Beam diameter 

was also varied between 0.63 5 and 0.07 cm. For each of the above cases, 

We 2 < 10" 2 

therefore, the momentum force of the electron beam will not be a primary 
cause of weld pool instability in either M551 or M553. 
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Appendix C 

This appendix contains detailed results of calculations for the M552 ex- 
periment which were performed during the Phase B period of this study. The 
first portion of the appendix addresses the capillary flow analysis, while the 
second segment discusses bubble dynamics associated with M552. 

C.l M552 CAPILLARY FLOW 

The effect of gap width and low gravity on capillary flow in the M552 
Exothermic Brazing Experiment has been studied. Time-to-spread, flow 
velocities and Reynolds number have been calculated for the configuration 
shown in Fig. C-l. The results are shown in Tables C-l, C-2 and C-3. 

The results indicate that: 

• Spreading time will be up to 50% shorter in gravity 
fields anticipated for Skylab than in ground tests. 

• Turbulent or oscillatory laminar flow will occur in 
Skylab processing in certain gaps where only laminar 
flow could occur in ground tests. 

The flow equation for an annulus is given by, 

u = (h o - h > < ai > 

where 

u = velocity 

h = axial position in annulus 
P = density 
R = sleeve radius 
/i = viscosity 

k = tube radius/ sleeve radius 

C-l 
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Table C-l 

M552 SPREADING TIMES 


Gravity 

(g/g E ) 

Gap Width 
(in.) 

Spread Time (Annulus) 
(jUsec) 

2 .'lO“ 6 

.005 

12.00 


.010 

6.00 


.020 

2.97 


.030* 

1.76 

1 

.005 

13.60 


.010 

7.47 


.020 

4.42 


.030 

3.60 


Assuming no variation in gap width. 


Table C-2 

M552 FLOW VELOCITIES 


Geometry 

Gravity 

(g/g E ) 

Gap 

(in.) 

Velocity (u) 
(cm/ sec) 

Annulus 

2 • 10 -6 

.005 

34 



.010 

67 



.020 

134 



.030* 

200 


1 

.005 

30 



.010 

55 



.020 

90 



.030* 

97 

Slot 

2 • 10' 6 

.005 

33 


\ 

.010 

66 



.020 

132 



.030* 

198 


1 

.005 

30 



.010 

53 



•020* 

81 



.030 

84 


* 

Assuming no variation in gap width. 
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Table C-3 

M552 REYNOLDS NUMBERS 


Geometry 

Gravity 

<g/g E ) 

Gap 

(in.) 

Re^ 

Flow 

Regime 

Annulus 

2-10' 6 

.005 

202 

Laminar 



.010 

794 

Oscillatory 



•° 20 * 

3200 

Turbulent 



.030 

7210 

Turbulent 


1 

.005 

182 

Laminar 



.010 

650 

Laminar 



.020 

2170 

Oscillatory 



.030* 

3490 

Turbulent 

Slot 

2*10' 6 

.005 

99 

Laminar 



.010 

395 

Laminar 



.020 

1580 

Laminar 



.030* 

3550 

Turbulent 


1 

.005 

90 

Laminar 



.010 

315 

Laminar 



.020 

970 

Laminar 



.030 

1500 

Laminar 


Assuming no variation in gap width. 
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h 


1 -k 4 + 1 - k 2 


1 -W 


fnk 


gravity 

2cr cos 8 
pgR(l-k) 


<y = surface tension 
9 = contact angle 


Equation (C.l) was developed assuming laminar, steady state, incompressible, 
isothermal flow with g opposing the flow direction. It should be a valid 
approximation for the M552 configuration and processing conditions. 


For very narrow gaps, Eq. (C.l) reduces to that for flow in a slot which 
is given by 


where 


u 


_ 2 

P g .P... 

12/xh 


(h o - h) 


D = slot width 


(C.2) 


This equation was developed for assumptions similar to those of Eq. (C.l) and 
will be just as rigorous except for geometrical effects. 

The spreading time, t, for flow in an annulus can be given by 


t = ft (h £n(l -h/h ) + h) (C.3) 

PgR 0 ° ° 

Again, this equation is valid for isothermal, laminar, incompressible flow 
with "g" opposing flow. Levich (Ref. C.l) states it will only be an approximation 

if t <^i . 

M 

The Reynolds number for flow in an annulus is given as (Ref. C-2), 
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or by (Ref. C-3), 


Re = 

1 11 


4R up 

Re = 

2 ii 


(G.4) 


(C.5) 


where R^ = hydraulic radius. 

The transition from laminar to turbulent flow occurs at Re^ ^ 2300. The 
transition from laminar viscous flow to oscillatory flow is given by Re^ =700. 
By making the transformation, k = 1 -e, and allowing e — *■ 0, it can be shown 
that Rej — " 11-^2 as k — *■ 1. Thus the transition points for Re^ are approx- 

imately 700 and 2300 also. The Reynolds number for flow in a slot is defined 
as 


Re = -£25 

» R 

' ■ 

and the laminar-turbulent transition is at Re ss 2000. 

s 

C.2 M552 BUBBLE DYNAMICS 

Many ground and KC-135 aircraft tests for M552 have exhibited a bubble 
or void space within the braze structure. Analysis of the movement of such 
a bubble or void in low gravity and the nonisothermal conditions of M552 has bee 
been initiated. 

Neglecting wall effects, the limiting velocity of a single bubble moving 
in a homogenous liquid under the combined influence of a gravity field and a 
temperature gradient is given by Ref. C-4, 
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u = 


3 u B 

4 rr 


. da dT , . 2 / . 

M g aA Tf d?- (Pg-P«> 8 a ( V M < 


where 


A = 2 + kVk 

* g 

B = 3 ii +2 i±. 

g £ 

k = thermal conductivity 
Ii = viscosity 
a = bubble radius 
c = surface tension 
T = temperature 
p = density 
g = gravity 


and subscripts £ and g refer to the liquid and gas phases, respectively. The 
temperature gradient affects bubble motion through tangential stresses at the 
bubble -liquid interface. These stresses are caused by surface tension gradients 
across the bubble. The surface tension mechanism is termed the Marangoni 
effect, which can be caused by either temperature or concentration gradients. 

The ratio of buoyant/Marangoni forces on the bubble is given by. 


y = 


P g ) ga (jU g + M £ ) A 

*g dT dz 


where y = buoyancy/surface tension forces. Values of y for various size 
bubbles and gravity levels are given below for M552 conditions using a 10°C/ 
cm temperature gradient. 
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Bubble 

Diameter 

(cm) 

Gravity 

(g/g E ) 

y 

io " 4 

1 

io 6 

io ' 4 

IO -6 

i 

io " 2 

1 

10 5 

io " 2 

10" 6 

io -1 


Thus, Marangoni forces will be important in Skylab, but not on ground tests. 
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Appendix D 

This appendix contains detailed results of calculations for the M553 
experiment which were performed during the Phase B portion of this study. 

The first section of this appendix contains analyses of the driving forces, 
magnitude and pattern of natural convective fluid flows which occur in the 
Sphere Forming Experiment. Following sections address such topics as: 
simple thermal analyses; vaporization; and molten pool splattering due to the 
impinging electron beam. 

D.l CONVECTION ANALYSIS 

As in the other M512 experiments, the fluid dynamics of the M553 
experiment is an important factor in determining the quality of the final product. 
Flow patterns in the molten material are important because all of the sample 
materials have low entropies of fusion (Ref. D-l). Thus, their solidification 
(microstructure) is controlled by the rate of heat transfer removal (Ref. D-2), 
which changes with the fluid flow (Refs. D-3 and D-4). The degree of flow will 
also determine the amount of mixing attained. If no or little flow were present, 
all heavier components would segregate to the lower portion of the spheres on 
earth, but to a negligible extent in the microgravity environment of space. 

Fluid flow can also affect the shape and release of the specimen while it is 
retained on the ceramic holder. 

Application of dimensional analysis (Refs. D-5 and D-6) to the governing 
equations for EB melting, coupled with ground-based and KC-135 experiments, 
should enable prediction of the extent of reduction or increase of motion in the 
molten metal and/or the change in flow pattern in electron beam melting in 
space. Possible physical forces which could induce fluid flow in the M553 
experiment, and their causes, include: 
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• Effective Gravity Force : Resultant force on weld specimen 
due to earth’s gravity and centrifugal and coriolis forces of 
orbiting spacecraft 

• Eorentz Force : Electromagnetic forces induced by passage 
of the electron beam current through the specimen 

• Elect restriction: Stresses induced when electrical permit- 
tivity changes with density 

• Magnetostriction : Stresses induced when permeability 
changes with density 

• Electrostatic Force: Caused by presence of excess elec- 

trical charge (due to beam current and/or thermionic 
emission) 

• Surface Tension : Tangential stresses at vapor-liquid 
interfaces can be induced if surface tension depends on 
temperature and/or concentration. Surface tension will 
also cause pressure gradients across curved interfaces. 

• Density Differences Accompanying Phase Changes 

• Beam Force : Impinging electrons give up their momentum 

• Thermal Expansion : Dilation and compression of fluids 
whose density changes appreciably with temperature can 
induce fluid flow 

• Vibration : Uncontrolled movement due to engine operation, 
astronaut motion, particle impacts, etc. 

• Centrifugal and Coriolis : Generated by disk rotation 

• Vapor Pressure : Evaporating molecules impart momentum 
which leads to normal stresses at vapor-liquid interface. 

• Inertia Forces : Tend to sustain induced motions 

• Viscous Forces : Tend to resist driving forces. 


D.1.1 Governing Equations 

The preceding forces, which could influence fluid flow and solidification, 
appear, explicitly in the following conservation equations which apply to form- 
ation of a molten pool by electron beam heating. 
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Continuity (Mass Balance): 


If + V • (pvj = 0 


where 

p = density 

V = velocity vector 
t = time 

V = "grad" or "del" operator 

Momentum: 


P + P ( V • V) V = -VP + p V 2 V + pg 

Inertia Viscous Gravity 

Force Force Force 

+ J x B + p E 

e 

Lorentz and Electro 
magnetostriction 
Forces 

where 


Electrostatic 

Force 


P = pressure 

p = viscosity 
1 

g = gravity vector 
.a 

J = electrical current density 
B = magnetic flux density 
P e = excess charge density 

E = electric field density 
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Ene rgy: 


pC y || + pC y (V • V) T = k V 2 x + + pQ 

+ - (P T/K) V • y 

Thermal Expansion 
Source 


where 



temperature 

heat capacity ^ 

thermal conductivity 
viscous dissipation function 
internal heat sources 
"conduction" current = o(E + V x B) 
electrical conductivity 


I (-§£) 


I 

p 3P 


) 


T 


Maxwell's Equations: 

VXE = - BB/at 

vxh = j + aD/at 

V • B = 0 

V • D = p 


where 


D = displacement current 
H = Magnetic flux 
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Constitutive Rel at ions ; 

D 

— I 

B 

* 

J 

P 

where 

e 1 = permittivity (dielectric constant) 
H ' = permeability 

Vapor -Liquid Boundary Conditions : 


= e ' E 

= p ' H 

= ct'(E + V x B) + p e V 

= P o [l - PAT + KApJ 


Force Balance, 



n. + 
1 


Beam and Vapor 
Pressure Forces 


as 

ax. 


+ s 


[*i 



n. 

l 


Surface Tension Forces 






+ 



n. 


where P^’ \ p vf ^ are, respectively, the pressures, 

viscosities, and velocity components in Phases 1 and 2, S is the surface 
tension, R^ and R^ are the principal radii of curvature of the surface, n. 
(i=l, 2, 3) are the components of the unit vector normal to the surface and 
directed into the interior of Phase 1, and summation over a repeated index 
(k=l, 2, 3) is assumed. 
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Energy Balance, 

ea <T 4 - T 4 ) - k n • VT = Q__ 
s oo s EB 

where 


e = emissivity 
a = Boltzmann's constant 
T g = surface temperature 

= environment temperature 

n = outward unit normal vector to free surface 
Q^B = heat flux from electron beam 


Continuity of Velocity, 


n x (V (i) - V (2) ) = 0 
n * V (1) = n • V (2) = 0 


Solid- Liquid Boundary Conditions ; 

Along stationary interfaces the no- slip condition (all components of velocity 
vanish) holds; whereas at moving interfaces (along melting or freezing fronts), 
a one -dimensional material balance yields 

u i = t 1 - P S /P L ) d6/dt 

where 

u^ = fluid velocity normal to the interface 

6 = position of the interface relative to the 
origin of the spatial coordinate system 

p g = solid density 

= liquid density 
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This is the source of the density-difference-accompanying-phase-change force 
for fluid flow. Also for a melting interface, an energy balance yields: 

-kV (T l - t s ) = p s X f 

where 

T = liquid temperature at the interface 
T_ = solid temperature at the interface 

D 

A, = heat of fusion 

and 



The preceding equations are based on the following assumptions: 

• No influence of external fields on physical properties 

• No coupling between constitutive flux relations (e.g., no Soret effect) 

• Single -component, Newtonian fluid 

• Constant physical properties 

• Bulk coefficient of viscosity vanishes. 


They are written in general vector form wherever possible, since they also 
apply directly to M551. It should be noted that only the differences in physical 
properties and geometry will differ in a dimensional analysis of M551 versus 
M553. 

D.1.2 Dimensional Analysis 

The controlling physical forces can be determined by nondimensioning 
each of the preceding equations and performing an order -of -magnitude com- 
parison on the various dimensionless groups which result (Refs.D-5 and D-6). 
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The key to successful analysis is in choosing the proper reference values; 
i.e., since no freestream velocity exists, which forces do we equate to esti- 
mate a "characteristic" or "typical" velocity. Oioosing the proper character- 
istic velocity is very important, since the reference time, temperature, etc. 
usually depend on this velocity. 

Previous examinations of film showing both KC-135 flight experiments 
and ground tests for the M553 experiment indicated fluid velocities exceeded 
300 cm/sec (p. B-5 of Ref. D-7). Continued M553 film analysis has since shown 
that flow velocities for nickel specimens on both KC-135 and ground tests have 
been approximately 20 cm/sec. Of the 14 possible forces affecting electron 
beam melting (Ref. D-8), only those shown in Table D-l yield characteristic 
velocities of this order of magnitude. These driving forces consist of couplings 
involving surface tension, gravity of Lorentz (electromagnetic) forces with 
inertia forces. Previous studies based on faster velocities had indicated coupling 
with viscous, rather than inertia, forces as controlling (Ref. D-7). This indi- 
cates the importance of choosing the correct characteristic velocity. Pre- 
liminary analysis also suggests that magnetostriction forces may also be im- 
portant, but gross uncertainty in electromagnetic property data for liquid metals 
precludes any decision at the present time. 

The dimensionless momentum equation which determines fluid flow in 
electron beam melting becomes (upon choosing the surface tension-inertia 
characteristic velocity). 


jT - pff + TT- p < V -V> V = -2jr“ VF+ V 2 V + ^(Pg)+-fJ MI (VxB)x B (D.l) 
N St dt 1N Oh Oh iN Oh 1N Oh 


where 


N 0 . = Stokes number = at / 0 L 

St o o 

. = duration of process/residence time 
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Table D-l 

PROBABLE CHARACTERISTIC VELOCITIES 
IN ELECTRON BEAM MELTING 



V elocity 

Controlling 
F orces 

Functional Form 

Value (cm/ sec) 
for Nickel in M553 

Inertia = Surface 
Tens ion 

(ST/pL) 1//2 

20 

■ 

Inertia = Lorentz 

(o' E 2 L/p) l/3 

1 

Inertia = Gravity 

(gL) l/2 

18 at 1 g E 
6 at 10 _4 g E 

Viscous = Lorentz 

(o' E 2 L 2 /p) l/2 

7 
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N 


Oh 


N 


Bo 


Ohnsorge number = fl/ ^ p^L S T q 
viscous force/surface tension force 

'y 

Bond number = Q g L /ST 


= gravity force/ surface tension 


N 


B 2 L 
o 


MI 


= Magnetic interaction number = -j — 7 g 7 p 
= magnetic force/surface tension 


o o 


Using physical property data for nickel and beam parameters for the 
M553 experiment, the equation reduces to 


10 3 p Jj +10 3 p (V -V) V = -10 3 VP + V 2 V + G(pg) + 10" 1 (V x B) xB 


where 

3 

G = 10 for earth gravity 
G = 10 * for expected Skylab gravity. 

This order-of-magnitude analysis indicates that surface tension driven con- 
vection will occur both in ground tests and for Skylab conditions because Nq^ 

« 1. This is confirmed by KC-135 M553 tests. Furthermore, gravity driven 
convection will exist on ground tests, but will be negligible in the reduced 
gravity of Skylab. Thus different forces will control convection on earth versus 
Skylab. The preceding analysis also indicates that electromagnetic or Lorentz 
forces will be negligible with regard to causing fluid motion. 

The surface tension driving force considered above is actually a surface 
tension gradient caused by radial and lateral temperature gradients. Ignoring 
convection, gradients of at least several hundred degrees Celsius per centi- 
meter have been predicted during melting (Ref. D-9). It can be shown by 
dimensional analysis (equate inertial and viscous terms) that the motion 
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caused by the initial temperature gradients will occur in less than 0.1 second 
and can persist for 60 seconds after removing the driving force. This means 
that there will be some fluid motion during solidification if the M553 specimens 
freeze after 30 to 40 seconds as predicted. 

Values of the pertinent dimensionless groups for the remaining M553 
materials are given in Table D-2. As can be seen, no significant changes 
are evident from the preceding conclusions for pure nickel. 

It should also be mentioned that the effects on fluid motion of surface 
free charges (excess electrostatic charge) and related electric and magnetic 
forces at the drop surfaces have been examined. The results of the examina- 
tion indicate that these surface charges are negligible driving forces for fluid 
motion. These results were also reviewed by Professor J. R. Melcher of MIT's 
Electrical Engineering Department, and he was in full agreement with these 
results (Ref. D-10). Furthermore, from the manner in which patches of surface 
contaminants moved about in the ground films, it is apparent that surface ten- 
sion driven flows due to concentration gradients are also important in this 
experiment. No reliable data exist, however, on what these impurities are, 
nor what the value of surface tension gradient with composition is. Thus the 
Marangoni effects are limited to thermal differences only in this study. 

D.1.3 Flow Patterns 

The aforementioned difference in controlling forces between terrestrial 
and space processing of the M553 experiment should result in altered flow 
patterns between these two conditions. From a knowledge of flow patterns 
generated by surface tension versus buoyancy controlled natural convection 
on Earth, the situation illustrated in Fig. D-l would no doubt exist during 
melting on Earth, while that shown in Fig. D-2 would prevail in space. 

The patterns shown in Fig. D-2 are based on the fact that surface tension 
is stronger in cooler regions; therefore the coolest portion of the specimen 
will exert the strongest "pull" on the surface. The buoyancy modification 
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Table D-2 

M553 DIMENSIONAL ANALYSIS 


Material 

N Oh 

Bo 

1 

* 

N Bo 

N t 

Bo 

N Oh 

N Oh 

N Oh 

Ni-Sn 

4.4 x 10" 4 

1.50 

2292 

3437 

.34 

Ni-Ag 

5.0 x 10' 4 

1.53 

1991 

3043 

.30 


-4 





Ni-Cu 

4.45 x 10 

1.54 

2107 

3241 

.32 


With Earth gravity (9.8 m/sec ) 
^With 10 4 Earth gravity 
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Specimen 



Fig. D-l - Convective Patterns in the M553 Specimen During 
Melting on Earth 



Fig. D-2 - Convective Patterns in the M553 Specimen During 
Melting in Space 
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exhibited in Fig.D-1 reflects the rise of warmer portions of fluid in a gravity 
field coupled with surface tension flows. 

D.2 M553 THERMAL ANALYSIS 

A critical question which must be answered for the M553 experiment 
is, "How long will it take for the released specimens to cool down to a "touch" 
temperature of 105°F?" Two separate models were formulated to study this 
problem. The results are shown in Table D-3. These results were based on 
the specimens being at the melting point {no superheat) upon release. 

Specimen conditions reported in Table D-3 include: free floating in 
vacuum; free floating and chamber backfilled with cabin atmosphere (5 psi 
mixture of 70% and 30% N^) two minutes after specimen release; and free- 
floating specimen and the chamber pumped with a steady stream of cabin at- 
mosphere at a flow rate of 5 cu ft/min. In the last case, the flow is considered 
to start two minutes after specimen release. 

Results shown in Table D-3 indicate that the chamber will be safe to open 
eight minutes after release of the last M553 specimen if a cabin atmosphere 
flow of 5 cu ft/min is established within the chamber two minutes after release 
of this last specimen. Otherwise, up to several hours may be needed for speci 
men cooldown to touch temperature. If the chamber is backfilled with 5 psi 
cabin atmosphere, cooldown time is roughly 35 to 40 minutes with 5 psi air 
versus several hours, unless the spheres remain attached to the chamber 
walls. It is unlikely, however, that all 11 of the released specimens will 
remain attached to a solid surface. 

The first model formulated was for a free-floating sphere being cooled 
by conduction and convection to an ambient atmosphere and/or being cooled by 
radiation to the chamber walls. Assuming that only sensible heat is being 
removed (no phase change) and that thermal diffusivity is infinitely large (an 
extremely good conductor), the energy balance becomes 


D- 14 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TH D306697 


Table D-3 

M553 COOLDOWN TIMES 


Thermal Condition 
of Specimen 

Cooldown 

Time* 

Chamber Wall 
Temperature ( F) 

Free floating in vacuum 

ESS 

80 



60 


E89 

-460 

Sticks to wall 

6 sec 

80 

Sticks to wall for 0.5 sec 

3.3 hr 

80 

Free floating and chamber 
backfilled with 5 psi 
mixture of 70% CL and 
30% N 2 

34 min 

80 

Free floating and 5 psi, 
70% O z - 30% N 2 cabin 

air flow at 5 cu ft/min 
through chamber 

7.7 min 

80 


For individual 0.25 in. molten nickel spheres at their melting 
point to reach 105°F. 
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where 


AT 4 4 

pV C Jt = -Aae (T 4 - T’ ) - Ah (T - T o ) (D.2) 


T = sphere temperature 
= wall temperature 

T q = chamber atmosphere temperature (% T ) 
a = Boltzmann's constant 
e = emissivity of sphere 
h = convective heat transfer coefficient 
A = surface area of specimen 
p = specimen density 
V = volume of specimen 
C = heat capacity of specimen 
t = time 


The coefficient h, for stationary cabin atmosphere, can be estimated by Ref. D-ll 


where 


k , Q , _ 1/4 1/3. 

= - ( 2 + Gr ' Pr ' ) 


k = thermal conductivity 
D = sphere diameter 
Gr = G rash of number 
Pr = Prandtl number 


For the condition of cabin air through the chamber, the heat transfer coefficient 
for the sphere (forced convection cooling of small spheres) can be obtained 
from a plot of Nusselt number versus Reynolds number and Prandtl number 
given in Ref. D-ll. 
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Equation (D.2) is valid for either a molten or solidified specimen. A 
computer program written in FORTRAN IV was written to solve Eq. (D.2) 
utilizing a fourth -order Runge-Kutta numerical algorithm. The output of this 
program is a temperature versus time profile. 

An estimate of the freezing time, t^, for this model, again assuming 
infinite thermal diffusivity is given by 


pRA 

3( ecI T 4 -h(T-T w )] 


(D.3) 


A correction time, t , for the effect of finite thermal conductivity is given by 

t = pCR 2 /k (D.4) 

Thus the time for a molten, free-floating specimen at its melting point to cool 
down to touch temperature is given by 


= t, + t 
f s 


+ t 


where 


t^. = total time to reach touch temperature 

t = time to remove sensible heat (from Eq. (D.2)) 
s 

The second model was developed for a specimen attached to the chamber 
wall being cooled by conduction to the wall and radiation. The chamber atmo- 
sphere was assumed to be a hard vacuum. Again, assuming that only sensible 
heat is being removed and that the specimen is a perfect conductor, the energy 
balance becomes, 

pVC 45" = -A' e <7 (T 4 - T 4 ) - A h (T - T ) (D.5) 

r 0t w w e w ' ' 
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where 


w 

A' 

h 


contact area between the specimen and chamber 
wall 

free surface area of specimen 

effective heat transfer coefficient to account for 
conduction losses 


Over the area of contact, A , the chamber wall can be assumed to be flat, 

w ’ 

therefore A^ and A' can be derived using simple geometry and are given 

b y 

A _ ?r 1//3 sin 2 e v 2//3 

w (1 - cos 9 ) 4//3 | 2 + cose j 2//3 


A' = 2 f l- l - cg i i-l A 
sin 2 0 1 W 


where 


0 = contact angle. 


The effective transfer coefficient is given by 

h = k /L 
e w 

where 


k = thermal conductivity of chamber wall 

L = some characteristic length for conduction 
in the wall. 


The characteristic length L can be estimated by 


L = Dk /2k 
w 
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The freezing time and correction time for finite conductivity can be 
estimated from Eqs. (D.3) and (D.4), respectively. 

Equation (D.5) can be solved using the algorithm developed for Eq. (D.2). 

D.3 VAPORIZATION IN M553 

A brief study of vaporization losses expected for nickel, aluminum, 
copper and silver under conditions expected during the melting portion of the 
M553 experiment indicates that aluminum will -evaporate 50 times faster than 
nickel, 25 times faster than silver and 10 times faster than copper. Thus, 
Ni-Af (5%) should not be considered as a material in the sphere forming ex- 
periment. The effect of concentration was included and the other sample 
materials were thus Ni-Ag(l%), Ni-Cu(30%), and pure nickel. 

The above data were obtained using Langmuir's formula 

” = P v tlsHjT k s/m 2 -sec 

where 

m = vaporization rate 
M = molecular weight 
P v = partial pressure 
T = surface temperature 
R = gas constant 

Surface temperatures were conservatively estimated to be 1725°K. 

D.4 MELT SPLATTERING INSTABILITY 

In both the M551 and M553 experiments, the heat released to the metal 
specimens by the impinging electron beam enables a molten pool to form. A 
question arises as to the stability of these molten liquids in low gravity. At 
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least one of the KC-135 M553 specimens was seen to break up violently into 
many smaller liquid spheroids upon complete melting. The electron beam 
was still hitting the specimen during break up. A similar instability might 
H develop in the molten puddle formed during the dwell mode of the M551 exper- 
iment, wherein the liquid might not adhere and separate from the solid disk 
(violent splattering). 

There are at least three different mechanisms which might explain the 
instability observed in the M553 KC-135 specimens. Upon melting, the specimen 
may have experienced violent degassing which could have led to droplet breakup. 
Another mechanism might be electrohydrodynamic instability caused by inter- 
actions between the electromagnetic forces for the electron beam and the fluid 
flow field set up in the molten metal by both thermal gradients, Lorentz and 
other forces. Lastly, the momentum force associated with the impinging 
electron beam might have set up unstable surface oscillations on the molten 
metal. 


The latter instability mechanism has been treated recently by Berghmans 
(Ref. D-12). He performed a theoretical study of fluid interface stability with 
special attention being given to the role of surface tension. A detailed examina- 
tion of this instability mechanism for both the M551 and M553 experiments is 
given in Appendix B where it is shown that the electron beam momentum force 
is not the cause of splattering in either of these experiments. 
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